Scale-free network topology and multifractality in weighted planar stochastic lattice 



O 

bj:): 

< 

o 



M. K. Hassan^ M. Z. Hassan^ and N. 1. Pavel^ 

1 Theoretical Physics Group, Department of Physics, University of Dhaka, Dhaka 1000, Bangladesh 

2 Institute of Computer Science, Bangladesh Atomic Energy Commission, Dhaka 1000, Bangladesh 

(Dated: August 31, 2010) 

We propose a weighted planar stochastic lattice (WPSL) formed by the random sequential parti- 
tion of a plane into contiguous and non-overlapping blocks and find that it evolves following several 



non-trivial conservation laws, namely 



-1 4/n- 



is independent of time V n, where Xi and yt 



are the length and width of the ith block. Its dual on the other hand, obtained by replacing each 
block with a node at its center and common border between blocks with an edge joining the two 
vertices, emerges as a network with a power-law degree distribution P{k) ~ where 7 = 5.66 
revealing scale-free coordination number disorder since P{k) also describes the fraction of blocks 
having k neighbours. To quantify the size disorder, we show that if the ith block is populated with 
Pi ~ Xi then its distribution in the WPSL exhibits multifractality. 
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I. INTRODUCTION 

Planar cellular structures formed by tessellation, tiling, 
or subdivision of plane into contiguous and nonoverlap- 
ping cells have always generated interest among scien- 
tists in general and physicist in particular because cel- 
lular structures are ubiquitous in nature. Examples in- 
clude acicular texture in martensite growth, tessellated 
pavement on ocean shores, agricultural land division ac- 
cording to ownership, grain texture in polycrystals, cell 
texture in biology, soap froths and so on For in- 

stance, Voronoi lattice formed by partitioning of a plane 
into convex polygons or Apollonian packing generated by 
tiling a plane into contiguous and non-overlapping disks 
has found widespread applications d, @ . There are some 
theoretical models, both random and deterministic, de- 
veloped either to directly mimic these structures or to 
serve as a tool on which one can study various physical 
problems However, cells in most of the existing struc- 
tures do not have sides which share with more than one 
side of other cells. In reality, the sides of a cell can share a 
part or a whole side of another cell. As a result, the num- 
ber of neighbours of a cell can be higher than the number 
of sides it has. Moreover, cellular structure may emerge 
through evolution where cells can be of different sizes and 
have different number of neighbours since nature favours 
these properties as a matter of rule rather than exception. 
A lattice with such properties can be of great interest as 
it can mimic disordered medium on which one can study 
percolation or random walk like problems. 

In this article, we propose a weighted planar stochastic 
lattice (WPSL) as a space-filling cellular structure where 
annealed coordination number disorder and size disorder 
are introduced in a natural way. The definition of the 
model is trivially simple. It starts with an initiator, say 
a square of unit area, and a generator that divides it 
randomly into four blocks. The generator thereafter is 
sequentially applied over and over again to only one of 
the available blocks picked preferentially with respect to 
their areas. It results in the partitioning of the square 




FIG. 1: A snapshot of the weighted planar stochastic lattice. 



into ever smaller mutually exclusive rectangular blocks. 
A snapshot of the WPSL at late stage (figure 1) pro- 
vides an awe-inspiring perspective on the emergence of 
an intriguing and rich pattern of blocks. We intend to in- 
vestigate its topological and geometrical properties in an 
attempt to find some order in this seemingly disordered 
lattice. If the blocks of the WPSL are regarded as isolated 
fragments then the model can also describe the fragmen- 
tation of a 2D object by random sequential nucleation of 
seeds from which two orthogonal cracks parallel to the 
sides of the parent object are grown until intercepted by 
existing cracks 0, S] • In reality, fragments produced in 
fracture of solids by propagation of interacting cracks is 
a formidable mathematical problem. The model in ques- 
tion can, however, be considered as the minimum model 
which should be capable of capturing the essential fea- 
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tures of the underlying mechanism. The WPSL can also 
describe the martensite formation as we find its defini- 
tion remarkably similar to the model proposed by Rao et 
al. which is also reflected in the similarity between the 
figure 1 and figure 2 of [H, Q . Yet another application, 
perhaps a little exotic, is the random search tree problem 
in computer science [lO[. 

Searching for an order in the disorder is always an at- 
tractive proposition in physics. To this end, we invoke 
the concept of complex network topology to quantify the 
coordination number disorder and the idea of multifrac- 
tality to quantify the size disorder of the blocks in WPSL. 
It is interesting to note that the dual of the WPSL (DW- 
PSL) obtained by replacing each block with a node at its 
center and common border between blocks with an edge 
joining the two corresponding vertices emerges as a net- 
work. The area of the respective blocks is assigned to the 
corresponding nodes to characterize them as their fitness 
parameter. Nodes in the DWPSL, therefore, are char- 
acterized by their respective fitness parameter and the 
corresponding degree k defined as the number of links a 
node has. For a decade, there has been a surge of inter- 
est in finding the degree distribution P{k) triggered by 
the work of A.-L. Barabasi and his co-workers who have 
revolutionized the notion of the network theory by rec- 
ognizing the fact that real networks are not static rather 
grow by addition of new nodes establishing links prefer- 
entially, known as the preferential attachment (PA) rule, 
to the nodes that are already well connected firl. Incor- 
porating both the ingredients, growth and the PA rule, 
Barabasi and Albert (BA) presented a simple theoreti- 
cal model and showed that such network self-organizes 
into a pow er-law degree distribution P{k) ~ k~^ with 
7 = 3 The phenomenal success of the BA model 

lies in the fact that it can capture, at least qualitatively, 
the key features of many real life networks. Interestingly, 
we find that the DWPSL has all the ingredients of the 
BA model and its degree distribution P{k) follows heavy- 
tailed power-law but with exponent 7 — 5.66 revealing 
that the coordination number of the WPSL is scale-free 
in character. 

In addition to characterizing the blocks of the WPSL 
by the coordination number fc, they can also be char- 
acterized by their respective length x and width y. We 
then find that the dynamics of the WPSL is governed by 
infinitely many conservation laws, namely the quantity 
^^x"~^y^^" ^ remains independent of time V n where 
blocks are labbeled by index i = 1, 2, N . For instance, 
total area is one obvious conserved quantity obtained by 
setting n = 2, sum of the cubic power of the length (or 

width) of all the existing blocks xf is a non-trivial 
conserved quantity obtained by setting n = 1 (or n = 4) . 
Interestingly, we find that when the «th block is pop- 
ulated with the fraction of the measure pi '-^ xf then 
the distribution of the population in the WPSL emerges 
as multifractal indicating further development towards 
gaining deeper insight into the complex nature of the 
WPSL we proposed. Multifractal analysis was initially 
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STEP 1 STEP 5 & 6 (dotted lines) 

FIG. 2: Schematic illustration of the first few steps of the 
algorithm. 

proposed to treat turbulence but later successfully ap- 
plied in a wide range of exciting field of research tl,3|. 
Recently though it has received a renewed interest as it 
has been found that the wild fluctuations of the wave 
functions in the vicinity of the Anderson and the quan- 
tum Hall transition can be best quantifled by using a 
multifractal analysis [13]. 

The organization of this paper is as follows. In section 
2, we give its exact algorithm of the model. In section 
3 various structural topological properties of the WPSL 
and its dual are discussed in order to quantify the an- 
nealed coordination number disorder. In section 4 we 
discuss the geometric properties of the WPSL in an at- 
tempt to quantify the annealed size disorder. Finally, 
section 5 gives a short summary of our results. 



II. ALGORITHM OF THE MODEL 

Perhaps an exact algorithm can provide a better de- 
scription of the model than the mere definition. In step 
one, the generator divides the initiator, say a square of 
unit area, randomly into four smaller blocks. We then 
label the four newly created blocks by their respective 
areas ai , 02 , 03 and 04 in a clockwise fashion starting 
from the upper left block (see figure 2). In each step 
thereafter only one block is picked preferentially with re- 
spect to their respective area (which we also refer as the 
fitness parameter) and then it is divided randomly into 
four blocks. In general, the jth step of the algorithm can 
be described as follows, (i) Subdivide the interval [0, 1] 
into (3j — 2) subintervals of size [0,ai], [ai,ai + 02], 

ai, I] each of which represents the blocks labelled 
by their areas ai, 02, 0(3^-2) respectively, (ii) Gen- 
erate a random number R from the interval [0, 1] and 
find which of the {3i — 2) sub-interval contains this R. 
The corresponding block it represents, say the pth block 
of area Op, is picked, (iii) Calculate the length Xp and 
the width Up of this block and keep note of the coor- 
dinate of the lower-left corner of the pth block, say it 
is {xiowtUIow)- (iv) Generate two random numbers xr 
and yR from [0, Xp] and [0, yp] respectively and hence the 
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FIG. 3: Shown is the growth of the number of blocks Nk{t) 
with exactly A; = 4, 5, 6, 7 neighbours as a function of time. 



point {xr + xiow-iUr + yiow) mimics a random point cho- 
sen in the block p. (v) Draw two perpendicular lines 
through the point {xji + xiow, Vr + yiow) parallel to the 
sides of the pth block in order to divide it into four smaller 
blocks. The label Qp is now redundant and hence it can be 
reused, (vi) Label the four newly created blocks accord- 
ing to their areas Up, a(3j_i), asj and 0(3^+1) respectively 
in a clockwise fashion starting from the upper left corner, 
(vii) Increase time by one unit and repeat the steps (i) - 
(vi) ad infinitum. 



III. TOPOLOGICAL PROPERTIES OF WPSL 

We first focus our analysis on the blocks of the WPSL 
and their coordination numbers. Note that for square 
lattice, the deterministic counterpart of the WPSL, the 
coordination number is a constant equal to 4. However, 
the coordination number in the WPSL is neither a con- 
stant nor has a typical mean value, rather the coordi- 
nation number that each block assumes in the WPSL is 
random. Moreover, it is allowed to evolve with time and 
hence the coordination number disorder in the WPSL 
can be regarded as of annealed type. Defining each step 
of the algorithm as one time unit and imposing periodic 
boundary condition in the simulation, we find that the 
number of blocks Nk{t) which have coordination num- 
ber k (or k neighbours) continue to grow linearly with 
time Nk{t) = nikt (see figure 3). On the other hand, 
the number of total blocks N{t) = '^/^Nkit) at time t 
in the lattice also grow linearly with time N{t) = 1 + 3t 
which we can write N{t) ~ 3t in the asymptotic regime. 
The ratio of the two quantities Pk{t) = Nk{t)/N{t) that 
describes the fraction of the total blocks which have co- 
ordination number k is pk{t) ~ mk/3. It implies that 
Pk (t) becomes a global property since we find it indepen- 
dent of time and size of the lattice. We now take WPSL 
of fixed size or time and look into its structural proper- 
ties. For instance, we want to find out what fraction of 
the total blocks of the WPSL of a given size has coor- 
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FIG. 4: Degree distribution P{k) for the DWPSL network 
where data points represent average of 50 independent real- 
izations. The line have slope exactly equal to 5.66 revealing 
power-law degree distribution with exponent 7 = 5.66. 
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FIG. 5: Cumulative degree distribution P{k' > k) is shown 
using using the same data as of figure 4. The dotted line with 
slope equal to 7 — 1 = 4.66 is drawn to guide our eyes. 



dination number k. For this we collect data for pf as a 
function of k and obtain the coordination number dis- 
tribution function pt{k) where subscript t indicates fixed 
time. Interestingly, the same WPSL can be interpreted 
as a network if the blocks of the lattice are regarded as 
nodes and the common borders between blocks as their 
links which is topologically identical to the dual of the 
WPSL (DWPSL). 

In fact, the data for the degree distribution P{k) of 
the DWPSL network is exactly the same as the coordi- 
nation number distribution function pt{k) of the WPSL 
i.e., pt{k) = P{k). The plot of lnP(A:) vs ln(A;) in fig- 
ure 4 using data obtained after ensemble average over 
50 independent realizations clearly suggests that fitting 
a straight line to data is possible. This implies that the 
degree distribution decays obeying power-law 



P{k) ~ k- 



(1) 



However, note that figure 4 has heavy or fat-tail, bench- 
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mark of the scale-free network, which represents highly 
connected hub nodes. The presence of messy tail-end 
complicates the process of fitting the data into power- 
law forms, estimating its exponent 7, and identifying the 
range over which power-law holds. One way of reduc- 
ing the noise at the tail-end of the degree distribution is 
to plot cumulative distribution P{k' > k). We therefore 
plot \n{P{k' > k)) vs In(fc) in figure 5 using the same 
data of figure 4 and find that the heavy tail smooths out 
naturally where no data is obscured. The straight line fit 
of figure 5 has a slope 7 — 1 = 4.66 which indicates that 
the degree distribution (figure 4) decays following power- 
law with exponent 7 = 5.66. We find it worthwhile to 
mention as a passing note that the mean coordination 
number kP{k) reaches asymptotically to a constant 
value equal to 5.333. 

We thus find that in the large-size limit the WPSL 
develops some order in the sense that its annealed co- 
ordination number disorder or the degree distribution 
of its dual is scale-free in character. This is in sharp 
contrast to the quenched coordination number disorder 
found in the Voronoi lattice where it is almost impossi- 
ble to find cells which have significantly higher or fewer 
neighbours than the mean value A; = 6 In fact, it 

has been shown that the degree distribution of the dual 
of the Voronoi lattice is Gaussian. The square lattice, on 
the other hand, is self-dual and its degree distribution is 
P{k) = 6{k — 4). Further, it is interesting to point out 
that the exponent 7 — 5.66 is significantly higher than 
usually found in most real-life network which is typically 
2 < 7 < 3. This suggests that in addition to the PA rule 
the network in question has to obey some constraints. 
For instance, nodes in the WPSL are spatially embedded 
in Euclidean space, links gained by the incoming nodes 
are constrained by the spatial location and the fitting 
parameter of the nodes. Owing to its unique dynamics 
this was not unexpected. Perhaps, it is noteworthy to 
mention that the the degree distribution of the electric 
power grid, whose nodes like WPSL are also embedded 
in the spatial position, is shown to exhibit power-law but 
with exponent 7powor = 4 [ll| . 

The power-law degree distribution P{k) has been 
found in many seemingly unrelated real life networks. 
It implies that there must exists some common under- 
lying mechanisms for which disparate systems behave in 
such a remarkably similar fashion f\3\ . Barabasi and Al- 
bert argued that the growth and the PA rule are the main 
essence behind the emergence of such power-law. Indeed, 
the DWPSL network too grows with time but in sharp 
contrast to the BA model where network grows by addi- 
tion of one single node with m edges per unit time, the 
DWPSL network grows by addition of a group of three 
nodes which are already linked by two edges. It also dif- 
fers in the way incoming nodes establish links with the 
existing nodes. 

To understand the growth mechanism of the DWPSL 
network let us look into the jth step of the algorithm. 
First a node, say it is labeled as a^,, is picked from the 



(3j — 2) nodes preferentially with respect to the fitness 
parameter of a node (i.e., according to their respective ar- 
eas). Secondly, connect the node Up with two new nodes 
(3j — 1) and (3j -I- 1) in order to establish their links with 
the existing network. At the same time, at least two or 
more links of Op with other nodes are removed (though 
the exact number depends on the number of neighbours 
Up already has) in favour of linking them among the three 
incoming nodes in a self-organised fashion. In the pro- 
cess, the degree kjp of the node Up will either decrease 
(may turn into a node with marginally low degree in the 
case it is highly connected node) or at best remain the 
same but will never increase. It, therefore, may appear 
that the PA rule is not followed here. A closer look into 
the dynamics, however, reveals otherwise. It is interest- 
ing to note that an existing nodes during the process gain 
links only if one of its neighbour is picked not itself. It 
implies that the higher the links (or degree) a node has, 
the higher its chance of gaining more links since they 
can be reached in a larger number of ways. It essen- 
tially embodies the intuitive idea of PA rule. Therefore, 
the DWPSL network can be seen to follow preferential 
attachment rule but in disguise. 

IV. GEOMETRIC PROPERTIES OF WPSL 

We again focus on the blocks of the WPSL but this 
time we characterize them by their length x and width 
y instead of the number of neighbours they have. Then 
the evolution of their distribution function C{xj_u[t) can 
be described by the following kinetic equation [3, Q 

dC{x,y;t) 

— = _xyC{x,y;t)+ (2) 

4 / / C{xi,yi;t)dxidyi. 

Jx Jy 

Incorporating the 2-tuple Millen transform 

/"OO /"OO 

M{m,n;t)^ / x'''-^y"~'^C{x,y;t)dxdy, (3) 
Jo Jo 

in the kinetic equation yields 

dM(m,n;t) / 4 \ , ^ , ^ 

Iterating it to get all the derivatives of M{m, n; t) and 
then substituting them into the Taylor series expansion 
of M{m, n; t) about t = one can write its solution in 
terms of generalized hypergeometric function [l7j 

M{'m,'n]t) — 2F2{^o,+ TO,-',Tn,n]—tj, (5) 

where M{m, n; t) = M{n, m; t) for symmetry reason and 
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FIG. 6: The plots of x'^-^y'^^^^'^ vs iV forn = 3, 4, 5 are 
drawn using data collected from one realization. 
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FIG. 7: The plots of a;f vs for four different realizations 
shows that the numerical value is different at every indepen- 
dent realization. 



One can immediately see that (i) M(l, l;t) = 1 + 3t is 
the total number of blocks N{t) and (ii) M(2, 2;t) = 1 
is the total area of all the blocks which is obviously a 
conserved quantity. The behaviour of M{m,n;t) in the 
long time limit is 



M{m,n;t) T 



(7) 



It implies that the system is in fact governed by several 

conservation laws namely, M(n, 4/n; t) = J^i' ^i~'^yr ^ 
are independent of time for all n which has been con- 
firmed by numerical simulation (see figure 6). 

We now focus on the distribution function n{x,t) — 
Jp°° C{x, y, t)dy that describes the concentration of blocks 
of length X at time t regardless of the size of their widths 
y. Then the gth moment of n{x,t) is defined as 



M,{t) 



x'^n{x, t)dx. 



(8) 



Appreciating the fact that Mq{t) — M{q + 1,1; and 
using equation ([7]) one can immediately write that 



Mq{t) - t 



{VV+16-(g+2)}/2 



(9) 



Note that had we focus on n(y,t) instead of n{x,t) we 
would have exactly the same result since their qth mo- 
ments are identical, M{q + 1, 1; t) — M(l, 1; t), due to 
symmetry reason. We find that the quantity M^lt) and 



hence J2i is a conserved quantity. How- 

ever, yet another interesting fact is that although Af3(t) 
remains a constant against time in every independent re- 
alization, their exact numerical value fluctuates sample 
to sample (see figure 7). It clearly indicates lack of self- 
averaging or wild fluctuation. Nonetheless, during each 
realization we can use M3{t) as a measure to populate 
the ith block with the fraction of the total population 
Pi — xf/J2i^i- The corresponding "partition function" 
then is 



(10) 



whose solution can be written immediately from equation 
dSl) to give 



Zq{t) - t{VV + 16-(3q+2)}/2^ 



(11) 



We find it instructive to express Zq in terms of square 
root of the mean area 6{t) = ^M(2, 2; f)/M(l, l;t) - 
as it gives the weighted number of squares N{q,S) 
needed to cover the measure which scales as 



where the mass exponent 



T{q) = VV + 16- (3g + 2). 



(12) 



(13) 



The non-linear nature of T(q) suggests that an infinite 
hierarchy of exponents is required to specify how the mo- 
ments of the probabilities {p}s scales with S. Note that 
t(0) = 2 is the Hausdorff-Besicovitch (H-D) dimension of 
the WPSL since the bare number N{q = 0, J) ~ (5~'^(°) 
and r(l) — O follows from the conservation laws (or nor- 
malization of the probabilities = 1). 

We now perform the Legendre transformation of T{q) 
by using the Lipschitz-Holder exponent 



a 



drjq) 
dq 



(14) 



as an independent variable to obtain the new function 



/(a) = qa + riq). 



(15) 



Substituting it in equation ((T2|) we find that 

7V(g,(5) - (16) 
since it is the dominant value of the integral 



Nia,S) - / pia)daS-f^"^6'^'^, 



(17) 
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It is interesting to note that ai is related to the gener- 
alized dimension Dq, also related to the Renyi entropy 
Hq [p) — -zzj In Si Pi in the information theory, given by 



Dr. — lim 



1 lnE.P? 



q-l In 5 



1-9' 



(22) 



which is often used in the multifractal formalism as it 
can also provide insightful interpretation. For instance, 
-Do = ''"(0) is the dimension of the support, Di — a\ is 
the Renyi information dimension and is known as the 
correlation dimension Il9l [201 • 



FIG. 8: The /(a) spectrum. 



obtained by the extremal conditions. We can infer from 
this that the number of squares dN{a, 5) needed to cover 
the PSL with a in the range a to a + da scales as 



(18) 



where p{a)da is the number of times the subdivided re- 
gions indexed by a [isj . It implies that a spectrum of 
spatially intertwined fractal dimensions 



/(«('?)) 



16 



16 



(19) 



are needed to characterize the measure which is always 
concave in character (see figure 8). It implies that the size 
disorder of the blocks are multifractal in character since 
the measure {pa] is related to size of the blocks. That 
is, the distribution of {pa} in WPSL can be subdivided 
into a union of fractal subsets each with fractal dimension 
f{a) < 2 in which the measure scales as 5". Note that 
/(a) is always concave in character (see figure 8) with a 
single maximum at (7 = correspond to the dimension of 
the WPSL with empty blocks. 

On the other hand, we find that the entropy S{5) — 
— Pi hipi associated with the partition of the measure 
on the support (WPSL) by using the relation X^iP? ^ 
in the definition of S{5). Then a few steps of 
algebraic manipulation reveals that S{5) exhibits scaling 



5(5) =ln5-"i 
ai = I obtaine( 
aq = - dT{q)/dq\^ . 



where the exponent ai = | obtained from 



(20) 



(21) 



V. SUMMARY 

We have proposed and studied a weighted planar 
stochastic lattice (WPSL) which has annealed coordi- 
nation number and the block size disorder. We have 
shown that the coordination number disorder is scale- 
free in character since the degree distribution of its dual 
(DWPSL), which is topologically identical to the network 
obtained by considering blocks of the WPSL as nodes 
and the common border between blocks as links, exhibits 
power-law. However, the novelty of this network is that 
it grows by addition of a group of already linked nodes 
which then establish links with the existing nodes follow- 
ing PA rule though in disguise in the sense that existing 
nodes gain links only if one of their neighbour is picked 
not itself. Besides, we have shown that if the blocks of 
the WPSL are characterized by their respective length 

X and width y then we find X^f^ ^^"^2/^^" ^ remains a 
constant regardless of the size of the lattice. However, 
the numerical values of the conserved quantities except 
n = 2 varies from sample to sample revealing absence 
of self-averaging - an indication of wild fluctuation. We 
have shown that if the blocks are occupied with a fraction 
of the measure equal to cubic power of their respective 
length or width then its distribution on the WPSL is mul- 
tifractal nature. Such multifractal lattice with scale-free 
coordination disorder can be of great interest as it has 
the potential to mimic disordered medium on which one 
can study various physical phenomena like percolation 
and random walk problems etc. 

NIP acknowledges support from the Bose Centre for 
Advanced Study and Research in Natural Sciences. 
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